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. ABSTRACT 

[•-[- ^ ■ Context. Recent studies have shown that stellar chromospheric activity, and its effect on convective energy transport in the envelope, 

I is most likely the cause of significant radius and temperature discrepancies between theoretical evolution models and observations. 

■ Accurate mass, radius, and abundance determinations from solar-type binaries exhibiting various levels of activity are needed for a 

■ better insight into the structure and evolution of these stars. 
Aims. We aim to determine absolute dimensions and abundances for the solar-type detached eclipsing binary EF Aqr, and to perform 
a detailed comparison with results from recent stellar evolutionary models. 

Methods, uvby light curves and uvbyP standard photometry were obtained with the Stromgren Automatic Telescope. The broadening 

■ function formalism was applied on spectra observed with HERMES at the Mercator telescope in La Palma, to obtain radial velocity 
curves. State-of-the-art methods were applied for the photometric and spectroscopic analyses. 

Results. Masses and radii with a precision of 0.6% and 1.0% respectively have been established for both components of EF Aqr The 
active 0.956 Mq secondary shows star spots and strong Ca 11 H and K emission lines. The 1.224 Mq primary shows signs of activity as 
^ ■ well, but at a lower level. An [Fe/H] abundance of 0.00+0. 10 is derived with similar abundances for Si, Ca, Sc, Ti, V, Cr, Co, and Ni. 

■ Solar calibrated evolutionary models such as Yonsei-Yale, Victoria-Regina and BaSTI isochrones and evolutionary tracks are unable 
' to reproduce EF Aqr, especially for the secondary, which is 9% larger and 400 K cooler than predicted. Models adopting significantly 
, lower mixing length parameters I /Hp remove these discrepancies, as seen in other solar type binaries. For the observed metallicity, 

Granada models with a mixing length of l/Hp = 1.30 (primary) and 1.05 (secondary) reproduce both components at a common age 
of 1.5±0.6 Gyr. 

Conclusions. Observations of EF Aqr suggests that magnetic activity, and its effect on envelope convection, is likely to be the cause 
of discrepancies in both radius and temperature, which can be removed by adjusting the mixing length parameter of the models 
downwards. 

vn 
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[ 1 . Introduction an overview. 'Clau sen et aH (l99^ and lTorres et al.1 (12006*) have 

I ■ shown that by reducing the mixing length parameter when calcu- 

^ . As recent studies suggest, current stellar evolutionary models lating the envelope c onvection, the model fits can be improved. 
_ ^ ; scaled to the Sun, are unable to predict the temperature and radii |ciausen et al.l (&) found that solar-type eclipsing binaries can 
of many active close bmary components with masses compa- (^jtjj exception) be divided in to two groups: slowly ro- 
_ rable to, or lower than the Sun. These models predict too high t^jj^g j^^^tj^g binaries, which can be fitted reasonably well by 
; temperatures and radii that are up to 10% lower than observed, j^^dels with a solar calibrated mixing length, and rapidly rotat- 
while the observations have errors as small as about 1 %, see e.g. j„g systems, which exhibit intrinsic variations and show signs 
[ Clausen et al.| (l2009l hereafter CBC09) and references therein. increased magnetic activity, which cannot be reproduced by 
Dynamo-generated magnetic fields can suppres convection and solar-calibrated models. These authors suggested that the evo- 
produce starspots, and are suggested as a cause for the discrep- j^jj^jj ^^e mixing length parameter with mass and activity is 
ancies in radii and temperature, see e.g. Ribas et al. (2008) for influenced by two different effects. There is a slight decrease in 

mixing length with increasing temperature and mass for inactive 

* Based on observations carried out at the Stromgren Automatic m^in sequence stars, and a strong decrease in the mixing length 
Telescope (SAT) at ESO, La Silla, and the Mercator Telescope, operated parameter for active stars compared to their inactive counterparts 
on the island of La Palma by the Flemish Community, at the Spanish c ■ ■^ n u - jj.ji 

^, ■ , , ,,,,,, r , T ■ , . r ■ of similar mass. However, more binaries are needed to develop 

Observatorio del Roque de los Muchachos of the Instituto de Astrotsica ,., . „ , ,. ^ ■ ■ , , -ii-i 

de Canarias ^ calibration for the evolution of mixing length with physical 

** T ui v/r u. fuc \ \- 1 1 ui ■ 1 f ■ f parameters and level of activity. 

Table X (light curves ot EF Aqr) is only available in electronic form ^ -' 

at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79. 128.5) or This paper is part of our ongoing study of several binary sys- 

via jhttp://cdsweb.u-s trasbg.fr/cgi-bin/qcat ?J/A-i-A/, tems with solar-type components, with either identical or differ- 
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ent levels of chromospheric activity (iClausen et al ] 120011 here- 
after CHOOl). In this paper the results of a complete analysis of 
the GO system EF Aqr (which has two components of different 
and relatively high activity) are presented. 



2. EF Aqr 

EF Aqr (my = 9.88, spectral type GO, 2'!85, also known as 
BD-7°5908, BV207, and as HD217512) is a solar-type eclips- 
ing binary that has not been studied in detai l before. It was firs t 
discovered at the Bamberg observatory bv IStrohmeiei ( 19581). 
In 1959 several times of minima were published by Nikulinal 
( 1959). Their studies do not agree well with the times of minima 
obtained for this paper, and were not used . In 2006 new ele- 
ments were published bv lOtero et al.l (l2006l) based on a combi- 
nation of data found in different catalogues. Both the published 
epoch and period differ from our results presented in Table [1] 
and Table[8] most likely because of the relatively low quality of 
the data Otero and collaborators were able to collect. In 2007, 
EF Aqr was included in the vari able star one sho t project, and 
one HARPS spectrum was taken dPall et alJl2007l) . It was anal- 
ysed as a single-lined spectrum, resulting in an effective temper- 
ature of 6110+187 K, a log § of 4.36 + 0.59 and a metallicity of 
[Fe/H]= 0.24 + 0.18. 



Table 1. Times of primary (P) and secondary (S) minima of EF 
Aqr determined from the uvby observations. 



3. Photometry 

Below, we prese nt the new photometric material for EF Aqr and 
refer to ICHOOll for more details on the observation and reduc- 
tion procedures, and determination of the times of minima. 



3.1. Light curves 

The differential uvby light curves of EF Aqr where obtained at 
the Stromgren Automatic Telescope (SAT) at ESO, La Silla and 
its six-channel uvby/3 photometer, on 80 nights between October 
2001 and December 2008 (JD 2452185 to 2454821). The ob- 
servations were made through an 18 arcsec diameter circular 
diaphragm. The light curves contain 764 points per band, and 
most orbital phases are covered at least twice. HD 217376, HD 
217877, and HD 218730 - all within a few degrees of EF Aqr 
on the sky - were used as comparison stars and were all found 
to be constant within a few mmag; see Table |2] The light curves 
are calculated relative to HD 217877, but all comparison star ob- 
servations were used, shifting them first to the same light level. 
The average accuracy per light curve point is 4-5 mmag (vby) 
and 5-6 mmag (m). The light curves (Table X) are available in 
electronic form at CDS. 

As seen in Fig.[Tl EF Aqr is a slightly active, detached system 
with quite different eclipse depths of about 0.7 mag and 0.2 mag 
(in y). The secondary eclipse occurs at phase 0.5 and the duration 
of the two eclipses is identical, meaning that the orbit is circular. 

3.2. Times of minima and ephemeris 

Four of each of the primary and secondary minima have been 
determined from the uvby light curves, and are listed in Table [1] 
By adopting a weighted least-squares fit to the times of pri- 
mary and secondary minima (i.e. assuming a circular orbit as 
argued for above), we derive the linear ephemeris given as 



HJD 
-2400000 


rms 


Type 


o-e 


52225.55854 


0.00014 


P 


0.00014 


52556.57275 


0.00010 


P 


-0.00001 


52576.54774 


0.00011 


P 


-0.00003 


52764.88340 


0.00100 


P 


-0.00012 


52563.70652 


0.00026 


S 


-0.00017 


52794.84594 


0.00014 


s 


-0.00009 


52917.54967 


0.00029 


s 


0.00004 


53676.59990 


0.00019 


s 


0.00011 



Mini = 2452556.57276(6)+2'!85357206(44)x£ 



(1) 



Notes. 0-C values are calculated for the ephemeris given in Eq. [T] 
assuming a circular orbit. 



3.3. Standard photometry 

Standard uvby/3 indices for EF Aqr and the c omparison stars, 
observed and derived as described by ICHOOlL are presented in 
Table|2] The indices are based on many observations and have a 
high precision. For comparison we included published photome- 
try from other sources as well. In general, the agreement is good, 
but in some cases differences are larger than the quoted errors. 
We have used only the new results from our own work for the 
analysis of EF Aqr presented here. 



3.4. Photometric elements 

Since EF Aqr is a well-detached system, JKTEBOP 
dSouthworth et alJ l2004allbl) has been used to determine the 
photometric elements of the uvby light curv es. JKTEBOP is 
based on the Nelson-Davis-Etzel binary model ( Nel son & DavisI 
1972; Etzel 1981; Martynov 1973), which represents deformed 
stars as biaxial ellipsoids and applies a simple bolometric 
reflection model. More details on the code and the general 
approach used here is explained in IClausen et a D (I2OO8, here- 
after CTB08). Hereafter the following symbols are used (with 
subscript s and p referring to respectively the secondary and the 
primary component of the system): / is the orbital inclination; 
e the eccentricity of the orbit; co the longitude of periastron; r 
the stellar radius in units of the semi-major axis; k - rjrp; u 
is the linear limb darkening coefficient; y the gravity darkening 
coeflicient; J the central surface brightness; L is the luminosity 
and Teff the effective temperature. 

The mass ratio between the components was kept at the spec- 
troscopic value (q = 0.7 61); se e section 14.21 Linear limb dark- 
ening coefficients from IClaretl (j2P00) were adopted. Extensive 
tests w ith nonlinear limb d arkening laws and Van Hamme coefli- 
cients dvan Hammelll993h were conducted, but showed little de- 
viation from the results presented here. The linear limb darken- 
ing law with Claret coeflicients was adopted because it resulted 
in the most consistent resul ts between the f our ba nds. Gravity 
darkening coeflicients from iBloemen et al. I (|20TTI) for convec- 
tive atmospheres were adopted. A circular orbit is assumed, fix- 
ing e cos w = e sin w = 0. We Justify this assumption in section 

The results for EF Aqr are presented in Table [3] The results 
from the four bands agree well. The errors quoted in this table 
are formal errors determined from the iterative least-squares so- 
lution procedure. Changing to VanHamme limb darkening coef- 
ficients, or adopting a nonlinear limb darkening law (quadratic. 
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Fig. 1. y light curve and b - y and u - b colour curves (instrumental system) for EF Aqr. 
Table 2. Photometric data for EF Aqr and the comparison stars. 



Object 


Spec. Type 


Ref. 


V 


cr 


b - y 


cr 


nil 


cr 




cr 


N(uvby) 




(T 




EF Aqr 


GO" 


Vll 


9.885 


9 


0.382 


5 


0.191 


8 


0.374 


10 


280 


2.621 


9 


8 






Vll 


10.092 


9 


0.356 


4 


0.178 


6 


0.386 


11 


22 








HD2 17376 


F3 IV" 


Vll 


6.807 


5 


0.288 


4 


0.149 


7 


0.478 


6 


192 


2.648 


7 


4 






083 


6.805 


6 


0.292 


4 


0.141 


6 


0.472 


6 


1 












094 


6.810 


4 


0.289 


3 


0.141 


4 


0.480 


5 


1 


2.645 


6 


1 


HD2 17877 


G2 V" 


Vll 


6.704 


4 


0.371 


4 


0.174 


7 


0.346 


7 


215 


2.608 


8 


4 


=HR8772 




P69 


6.65 




0.364 




0.181 




0.319 




2 












G076 


6.702 


6 


0.367 


4 


0.175 


6 


0.350 


22 


3 












G077 




















2.601 


8 


5 






AJ85 




















2.602 


2 


11 


HD2 18730 


Gl/2 V 


Vll 


7.333 


5 


0.385 


5 


0.198 


7 


0.361 


8 


180 


2.608 


7 


4 






083 


7.334 


6 


0.376 


4 


0.207 


6 


0.346 


6 


1 












094 


7.331 


4 


0.384 


3 


0.195 


4 


0.357 


5 


1 









Notes. For EF Aqr, the uvbyfi information by VI 1 is the mean value near phases 0.25 and 0.75 (first line), and at total secondary eclipse (second 
line). 

N is the total number of observ ations used to form th e mean values, and cr is the rms error (per observation) in mmag. 
Henry Draper Catalogue; "'' iHouk & Swifll ( ll999l) 

References. AJ85 = [Andersen & Jensenl ( 119851) : G076 = lGr0nbech & OlsenI ( fT97l : G077 = lGr0nbech & OlsenI ( Il977h : 083 = l01serj ( fl983l) : 
094 = (Mvfcv): l01sei] Jl994h . (BY. Olsen unpublished: P69 ^ iPerrvId 19691) : Vll = this paper 



logarithmic and square root) changes the resulting inclinations 
by about 0.1°, and the radii by 0.005. 

The limb darkening coefficients are based on the effective 
temperature of both components. The effective temperature was 
calculated from a weighted mean of three colour-metallicity- 
temperature calibrations, the spectroscopic temperature (see 
Section|5j, and for the secondary the flux scale method of Popper 
(119801) is included as well (see Section |6]l. To obtain the dered- 



dened colour indices, a reddening of £'(/7-3') = 0.0 18+0.01 1 wa s 
determined with the intrinsic colour cahbration o f l01senl(ll988h . 

The li ghtcti rves show clear indications for star spots. 
PHOEBE (Prsa 2006), a binary modelling code built on the 
Wilson-Devinney code (Wilson 2008), was used to analyse sep- 
arate sections of the light curves while including star spots in the 
model. Sections of different observation periods ranging from 
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Table 3. Photometric solutions for EF Aqr from the JKTEBOP 
code. 





u 


V 


h 


>' 


in 


88.40 


88.46 


88.44 


88.49 




±7 


± 6 


±5 


±5 




0.1224 


0.1213 


0.1218 


0.1212 




0.0875 


0.0870 


0.0873 


0.0876 


k 


0.7152 


0.7170 


0.7168 


0.7159 




±25 


± 19 


± 17 


± 15 


rp + r. 


0.2099 


0.2084 


0.2091 


0.2080 




+ 6 


± 5 


± 5 


± 5 


Up 


0.8107 


0.7673 


0.7369 


0.6553 




0.9223 


0.8796 


0.8323 


0.7478 


yp 


0.6527 


0.5648 


0.4465 


0.3791 




1.2437 


0.8670 


0.6904 


0.5580 


Js / Jp 


0.2438 


0.2851 


0.3729 


0.4251 




± 27 


+ 21 


± 19 


± 19 


/ ^p 


0.1184 


0.1410 


0.1834 


0.2092 


cr (mmag) 


10.9 


8.3 


7.5 


7.3 



Notes. Linear limb darkening coefficients from IClaretl l l2000h 
were adopted together with gravity darkening coefficients from 
iBloemen et alj li201 1) . A mass ratio of q = 0.761 from the spectroscopy 
was assumed, the orbit was assumed to be circular, and phase shift and 
magnitude normalisation were included as free parameters. 
The quoted errors are formal errors from the iterative least-squares so- 
lution procedure. 

Table 4. Adopted photometric elements for EF Aqr. 



i 


88.454 ± 0.074 




0.0 


0)' 


90.0 


rp 


0.1216 + 0.0012 




0.0871 + 0.0003 


rp + r. 


0.2098 ± 0.0013 


k 


0.7164 + 0.0063 





U 


V 


b 


y 


Js/Jp 


0.244 


0.285 


0.373 


0.425 


±9 


±9 


±12 


±13 


^s/ ^p 


0.1188 


0.1413 


0.1837 


0.2096 




+ 19 


±20 


±17 


±17 



Notes. adopted parameters 

The individual flux and luminosity ratios are based on the mean stellar 
and orbital parameters. 

five months to two years were tested, but PHOEBE was not able 
to fit the effect of the star spots. 

The adopted photometric elements listed in Table |4] are the 
weighted mean values of the JKTEBOP solutions in the four 
bands. Realistic errors were determined from Monte Carlo error 
simulations with 5000 iterations in each band and a comparison 
between the solutions in each band. The uncertainty of the limb 
darkening values based on the accuracy of the effective temper- 
ature was taken into account as well. The accurate light curves 
allow the relative radii to be determined with a precision higher 
than 1%. The obtained luminosity ratios correspond very well to 
those derived from the depths of secondary eclipse. 

4. Spectroscopy 

For the radial velocity and abundance determinations of 
EF Aqr, 17 high-resolution spectra were obtained with the 
HERMES spectrograph (R = 85000, 55 orders, 3770-9000 A, 



iRaskin et al.l (l201 ih X at the 1.2m Mercator telescope, Roque 
de los Muchachos Observatory, La Palma. The observations 
were obtained in October and November 2010 (JD 2455509- 
2455545). HERMES was used in high-resolution mode, and Th- 
Ar-Ne frames were obtained just before or just after the observa- 
tions to reduce the instrumental drift. The exposure times were 
adjusted to obtain a signal to noise ratio (S/N) around 100 in 
V. For the basic reduction of the spectra the HERMES pipeline 
version 3.0 was used together with dedicated IDE programmes 
to further clean and normalise the spectra. For the radial veloc- 
ity determinations only spectra from orders 13-48 (4016-6796 
A) were used, because the S/N in the other orders was too low. 

4. 1 . Radial velocities 

To obtain the radial velocities from the spectra, the broad- 
ening function (BE) formalism (Rucinski 1999, 2002, 200j) 
was used on 16 spectra. The last spectrum, which was ob- 
served just after secondary eclipse, was not included in the 
radial velocity determinations. Template spectra calculated for 
the primary and secondary components were used in the anal- 
ysis. The synthetic template spectra were computed with the 
bssynth tool (Bruntt, private communication), which applies 
the SYNTH software (Valenti & Piskunovlll996') and modified 
ATLAS9 models (Heiter et al. 2002). Regardless of the signifi- 
cant temperature difference between the components, both tem- 
plates give compatible results. But the radial velocities obtained 
with the primary template correspond slightly better to the fi- 
nal orbital solution. The broadening functions were calculated 
separately for the different orders, and then combined by using 
a weighted average based on the S/N level and the wavelength 
range of each order. 

The radial velocities were obtained from the broaden- 
ing functions by fi tting an analytical BE to the observed BE 
(iKaluznv et al.ll2()0 6). which was derived by assuming two per- 
fectly spherical rigid rotators. The assumptions made in the 
derivation of the analytical BE are valid for EF Aqr because both 
components show little or no deformation. Fitting Gaussians to 
the BF does not result in values significantly different from those 
presented in Tab. 5. 

The obtained radial velocities are given in Table |5] By ex- 
amining the distribution of the radial velocities calculated sepa- 
rately for each order, the uncertainties of the radial velocities are 
estimated at 0.35 km s"' for the primary and 1.31 km s"' for the 
secondary. 

4.2. Spectroscopic elements 

The spectroscopic orbits were calculated using the method of 
Lehman-F ilhes implemented in the SBOP programme (EtzeJ 
119811 12004 . which is a modified and extended version of an 
earlier code written bv i Wolfe et al.l (Il967h . The radial velocities 
are given a weight factor of 1 .0 except for those obtained from 
four spectra that had a significantly lower S/N value (taken at 
HJD=555 10.49097, 55545.38178, 55545.39625, 55545.40725). 
These radial velocities were given a weight factor of 0.5. To ob- 
tain accurate uncertainties, a wrapper with a Monte Carlo algo- 
rithm was written around the SBOP code. The input data were 
perturbed with the estimated uncertainties on the radial velocity 
measurements. The final uncertainties on the spectroscopic ele- 
ments are determined from 5000 iterations in the Monte Carlo 
code. 
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Table 5. The radial velocities of EF Aqr together with the resid- 
uals of the final spectroscopic orbit given in Table |6] Supscripts 
p and s denote determinations based on templates of the primary 
and the secondary. 



BJD 


Phase 


RV„ 


RV., 


(O-C)p 


(0-C),, 


-2400000 




km s"' 


km s"' 


km s"' 


km s"' 


55509.34736 


0.764 


79.69 


-114.33 


-0.08 


-0.03 


55509.38693 


0.778 


78.74 


-113.24 


-0.06 


-0.22 


55510.33598 


0.111 


-58.49 


67.06 


-0.45 


0.16 


55510.49097 


0.165 


-76.83 


91.89 


-0.28 


0.65 


55512.40508 


0.836 


67.96 


-98.88 


-0.20 


0.15 


55512.49166 


0.866 


58.41 


-86.06 


-0.25 


0.48 


55513.40673 


0.187 


-81.85 


98.23 


-0.14 


0.20 


55513.50138 


0.220 


-86.73 


104.43 


0.03 


-0.24 


55514.51857 


0.577 


35.21 


-55.25 


0.35 


0.00 


55516.43170 


0.247 


-88.10 


106.77 


0.13 


0.17 


55517.32517 


0.560 


27.23 


-44.85 


0.28 


-0.01 


55517.50485 


0.623 


55.03 


-81.42 


0.32 


-0.07 


55518.38558 


0.932 


30.76 


-50.13 


-0.22 


0.02 


55545.38178 


0.392 


-56.43 


64.66 


0.41 


-0.67 


55545.39625 


0.397 


-54.36 


61.59 


0.37 


-0.96 


55545.40725 


0.401 


-52.80 


59.65 


0.28 


-0.74 



Table 6. Spectroscopic orbital solution for EF Aqr obtained with 
the SBOP code. 



Parameter Primary Secondary 



p{dy 


2.85357206 


To" 


2452556.57276 


e" 











90 




0.761+0.003 


7 (km s"') 


-4.2341 + 0.0694 


-4.372 ± 0.0974 


K(km s-') 


84.2188 + 0.0882 


110.8169 ± 0.1238 


a sin i (Rq) 


3.3047 ± 0.0035 


4.3484 ± 0.0049 


M sin' / (Mo) 


1.244 ±0.008 


0.946 ± 0.006 


cr 


0.28 


0.39 



Notes. Assumed values. 

The quoted errors are formal errors from the iterative least-squares so- 
lution procedure. 



As in the determination of the photometric elements, a circu- 
lar orbit was assumed, fixing e = 0.0 and oj - 90°. The resulting 
spectroscopic elements are given in Table |6l The radial veloc- 
ities of the components were analysed independently (single- 
lined spectroscopic binary SB 1 ), but SB2 solutions lead to nearly 
identical results. 

In Fig |2] the spectroscopic orbital solution is plotted to- 
gether with the residuals of both components. As can be seen 
from the plot, there is a clear sine-like pattern with an ampli- 
tude of around 500 m s ' in the residuals of both components. 
Including the eccentricity and argument of periastron as free 
parameters improves the solution, but the reduced masses and 
semi-major amplitudes do not change. In this way an eccentric- 
ity of e = 0.0018 + 0.0022 and an argument of periastron of 
0) = 90.23° + 0.21 are found. However, the photometry (light 
curves and times of minima) limits the eccentricity to a max- 
imum of e = 0.00005 ± 0.00008, indicating that the orbit is 
practically circular. The periodic residuals in the radial veloc- 
ity curves can be explained by the activity of the system. The 
lig ht curves sho w clear signs of star spots, and previous research 
of iHuerta et alj (|2008) has shown that star spots can cause ra- 
dial velocity variations of the order of 0.5 km s This theory 




c I I 1 I I I I a I I I I I 1 I I I I I : 

0.0 0.2 0.4 0.6 0.8 1.0 



Phase 

Fig. 2. Spectroscopic orbital solution for EF Aqr (solid line: pri- 
mary, dashed line: secondary), and the observed radial velocities 
(filled circles: primary, open circles: secondary). Phase 0.0 cor- 
responds to the central primary eclipse. 

can not be proven because there are no accurate photometric ob- 
servations available that coincide with the spectroscopic obser- 
vations. Altogether, these small deviations in the radial velocity 
curves do not have a significant effect on the resulting spectro- 
scopic elements, and can be neglected. 

5. Chemical abundances 

To determine the chemical abundances of EF Aqr, we used all 
17 HERMES spectra. Because sufficient spectra are available, 
the DisEntangler package of Simon & Sturm (1994) was used 
to disentangle the spectra, so that the primary and secondary 
could be analysed separately. In this way a higher accuracy can 
be obtained than if composite spectra were used. The 35 orders 
used also for the radial velocity determinations, which had suf- 
ficient S/N and no influence of telluric lines, were disentangled 
individually for the abundance analysis. The DisEntrangler al- 
gorithm treats the spectra with an equal light ratio, resulting 
in component spectra that are determined independently of an 
additive constant that depends on the light ratio in each order. 
This constant is c alculated using MARCS model atmospheres 
dGustafsson et alj |2008). and implemented in the normalisation 
process using dedicated IDL programs. DisEntangler also as- 
sumes a constant light level, but since EF Aqr is constant to 
within 0.5 % outside eclipses, this is of no concern. A 40A region 
of the spectrum of the primary component, centred on 5585A is 
shown in Fig. [3] 

The ve r satile wavelength analysis (VWA) tool of 
iBruntt et"aD (l2002h was used to determine the abundance 
and temperature of bot h components of EF Aqr. VWA used 
the SYNTH software o f IValenti & Piskunovl (Il996h to generate 
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Fig. 3. 40 A region centred at 5585 A of the disentangled spectrum for the primary component of EF Aqr (black) and a matching 
synthetic spectrum calculated with VWA (blue). The lines indicated in green are used in the abundance determinations. 



synthetic spectra. Atmosphere models were i nterpolat e d from 
the MARCS model atmospheres (Gustafs son et alJ 



2008) 



which use the solar composition of iGrevesse et al.l (l2007h . 
The atomic line d ata were taken from the VALD database 
dKupka et al.ll 19991) . However, the log(.g f ) va l ues w ere adjusted 
so that every line measured bv lWallace et al.l ( 1998[) rep roduces 
the atmospheric abundances by Grevesse et alJ ( ^OOTh . For a 
detailed description of the VWA tool, the reader is refeiTed to 
[Brunttet al. (2004. 2008); Bruntt (2009); Bruntt et al. (2010a b). 

When comparing the observed orders to synthetic spectra, 
the rotational velocities of the two components are found to be 
vsinip -21 +4 km s" ' and v sin is = 18+4 km s" ' respectively, 
indicating that they are fast rotators. Owing to the higher degree 
of line blending caused by the fast rotation, only 468 (primary) 
and 396 (secondary) lines were selected. The surface gravity was 
fixed at the values determined from the photometric and spectro- 
scopic elements given in Section |6] Spectroscopy alone is usu- 
ally only capable to determine log ^ with ari accu racy of up to 0. 1 
dex (iBruntt et alJ2010atrvilenti & Fisched2005.) . but in a binary 
system accuracies as high as 0.01 dex can be obtained. The ef- 
fective temperature and micro-turbulence were tuned to remove 
differences between Fe i and Fe ii abundance, and any correla- 
tion between abundance and excitation potential or equivalent 
width. In this way effective temperatures of Teff.p - 6140 ± 50 
K and Tgff.s = 5200 + 100 K were found. iBruntt et al.1 (1201 Oah 
found a systematic offset of 40 K between spectroscopicaly ob- 
tained temperatures and temperatures obtained with direct meth- 
ods. After subtracting this offset spectroscopic temperatures of 
Jeff = 6100 + 50 K (primary) and Teff = 5160 + 100 K (sec- 
ondary) were obtained. 

The micro-turbulence is determined to be = 1.60 + 0.10 
km s ' and = 1.32 + 0.10 km s For the primary this 
agrees very well with the predicted micro-turbulence of - 
1 .65 ± 0.34 km s ' by the calibration of lEdvardsson et al.l (Il993h 
(which is an analytical fit to models based on 157 stars). For the 
secondary the predicted micro-turbulence is = 0.66 ± 0.38 
km s which is significantly lower than found in the abun- 
dance determination. However, the secondary's effective temper- 
ature is outside the range of the Edvardsson calibration, which 
is only valid in the temperature range T^ff - 5550 — > 6800 A'. 
IBruntt et"al] (1201 Oal) derived a micro-turbulence relation using 
23 solar-type stars. This calibration predicts ^ - 1 .27 ± 0.13 km 
s"' (primary) and ^ = 0.84 ± 0.1 1 km s"' (secondary). Both the 
predictions for the primary and the secondary are significantly 
lower than determined in the present work from observations. 



We cannot pinpoint an exact reason for the difference between 
the calibrations and the results presented here, but it has to be 
noted that the implementation of the microturbulence can de- 
pend on the model atmosphere used, which could explain the 
difference between our results and the results of Edvardsson. 

The abundances determined from all useful lines with an 
equivalent width between 10 and 90 mA are presented in Table 
|7] Obviously a very robust iron abundance of [Fe/H] - -0.002 
+ 0.090 dex is obtained with excellent agreement between the 
results of Fe i and Fe ii lines of both components. Changing 
the temperature of the primary by ±50 K and the secondary by 
+ 100 K results in a +0.04 dex change in Fe i abundance and 
+0.02 dex change in Fe n abundance. If the micro-turbulence is 
increased or decreased by +0.10 km s"', the neutral iron abun- 
dance changes by +0.07 dex, while the ionised iron abundance 
changes by +0.02 dex. The abundances of the other ions listed 
in Table |2] agree well with the results obtained for the iron abun- 
dance and a relative abundance identical to the solar value. The 
few outliers such as Cu i, C i, S i, and Ce n are based only on 
one line and can be neglected. The final metal abundance of EF 
Aqr is determined based on the abundances of all elements with 
eight or more accepted lines, and results in [M/H] = 0.00 ± 0. 10 
dex. 

In addition t o the spectroscopic abu ndance determinations, 
the calibration of iHolmberg et al.l (l2007h was used to determine 
the metalicity directly from the de-reddened uvby indices for the 
individual components. This results in metallicities of [Fe/H] = 
0.04 ± 0. 1 1 (primary) and [Fe/H] = -0.20 + 0.27 (secondary), re- 
spectively. The photometric metallicity obtained for the primary 
coiTesponds to the spectroscopic metallicity within errors. For 
the secondary the difference is larger, but the two methods still 
agree within the estimated errors. 



6. Absolute dimensions 

The absolute dimensions of EF Aqr are presented in Table [8] 
as calculated from the photometric and spectroscopic elements 
given in Tables |4] and |6] The accuracy for the masses is better 
then one percent while for the radii an accuracy of 1.2% was 
established. The V magnitudes and uvby indices were calculated 
from the combined magnitudes and indices of the system outside 
the eclipses (see Table |2]i. The V magnitude and uvby indices 
obtained for the primary component during secondary eclipse 
agree reasonably well with those obtained from the photomet- 
ric solution. The E{b - y) interstellar reddening was calculated 
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Table 7. The abundances ([El/H]) for the primary and secondary Table 8. Astrophysical data for EF Aqr 

















Systemic parameters 






Ion 


Primary 
[El/H] rms 


N" 


Secondary 
[El/H] rms 


N" 


P(d) 
To (HJD) 




2.85357206 
52556.5728 


Ci 


-0.24 




1 








e 




0.0 


Nai 


0.08 




1 


0.09 




1 


CO 




90.0 


Mgi 


0.01 




1 








a(Ro) 




10.9984 ± 0.0207 


All 


0.17 




1 








vo (km s"') 




-4.0660 + 0.1264 


Si I 


0.06 


0.08 


10 


0.10 


0.12 


3 


q 




0.761 + 


0.003 


Si II 


0.07 


0.09 


2 








i n 




88.454 = 


t 0.076 


Si 


0.35 




1 








E(b-y) 




0.018 + 0.011 


Cai 


0.01 


0.09 


5 


0.00 


0.19 


6 


[Fe/H] (dex) 




0.00 ±0.10 


Sc II 


-0.09 


0.11 


2 


0.12 


0.17 


2 


rf(pc) 




172.0 


±3.8 


Til 


-0.04 


0.10 


2 


-0.09 


0.15 


3 










Tin 


-0.06 


0.09 


3 


-0.03 




1 


Component parameters 






Vi 


-0.15 


0.11 


2 


0.24 


0.11 


8 




Primary 


Secondary 


Cri 


-0.09 


0.11 


9 


0.10 


0.11 


4 


M(Mo) 


1.244 d 


: 0.008 


0.946 ± 0.006 


Cm 


-0.02 


0.08 


4 










1.338 d 


: 0.012 


0.956 ± 0.012 


Mm 


-0.06 


0.12 


4 


0.18 


0.26 


4 


a: (km s"') 


84.180 d 


: 0.185 


110.686 ±0.302 


Fe I 


0.00 


0.08 


136 


0.00 


0.10 


88 


V sin I (km s"') 




21 ±4 


18 ±4 


Fe n 


0.00 


0.09 


11 


0.00 




1 


f (km s"') 


1.60 


± 0.10 


1.32 ± 0.10 


Co I 


0.13 


0.12 


2 


-0.06 


0.17 


2 


logg (cgs) 


4.280 d 


: 0.007 


4.453 ± 0.011 


Nil 


-0.06 


0.08 


18 


0.12 


0.15 


10 


Vo 


10.015 d 


: 0.022 


11.711 ± 0.024 


Cut 


-0.42 




1 


-0.62 




1 


(b-y) 


0.338 d 


: 0.010 


0.484 ± 0.016 


Zn I 


-0.19 




1 








m\ 


0.183 d 


: 0.010 


0.332 ± 0.021 


Yii 


-0.06 




1 


0.41 




1 




0.383 d 


: 0.015 


0.292 ± 0.029 


Cen 


-0.59 




1 








7^eff(K) 


6150 ± 65 


5185 ± 110 
















%ol 


3.847 d 


: 0.029 


5.319 ± 0.051 


Notes. Number of lines used 


per ion, 


the rms error is computed if 


log LI Lq 


0.360 d 


: 0.011 


-0.227 ± 0.020 


two or more lines are available. 










BC 


-0.021 d 


: 0.013 


-0.192 ±0.035 
















Mv 


3.868 d 


: 0.030 


5.511 ±0.0523 



using the calibration o f lOlsenI ([T988I) . using the ubvyP standard 
photometry for combined light outside eclipses. 

The temperature of both components was determined 
using the three different colour-metallicity-temperature cali- 
brations of lAlonso et al.1 (Il996i) . iHolmberg et alj (l2007l) and 
ICasagra nde et alJ ( 20101) . the spectroscopic derivation of iPoppeJ 
( 1980), and for the secondary also the flux scale method of this 
author. The results for each method are shown in Table |9l The 
uncertainties include the error on the colour indices, the metallic- 
ity, reddening and the uncertainty of the calibrations themselves. 
The obtained temperatures are mutually consistent within the er- 
rors of the calibrations, and as final temperatures the average val- 
ues over the different methods are adopted: T^jj p - 6150 + 65 
K and r,//,, = 5185 + 110 K. 

The projected rotational velocities from the broadening func- 
tions vsin/p = 23 km s"' and vsin/, = 17 km s"' agree 
with those from the abundance analysis in the previous section. 
The projected rotational velocities agree within the errors with 
the calculated velocities assuming a synchronous circular orbit: 
vsinip = 23.7 km s"' and vsin/j - 16.9 km s '. Furthermore, 
the age (1.5 ± 0.6 Gyr) determined in Section [8] is clearly older 
than the circula risatio n and synchronisation time scale as deter- 
mined by Za hr3(ll977h:T„>^ = 0.28 Gyr and t,v„ = 0.09 Gyr, and 
by Claret (iTorres et al.l2009l and references herein): Tc„c = 0.42 
Gyr, Tjy„(prim) = 0.21 and Tjv«(sec) - 0.20. 

The distance to the sy stem was calculated using the classical 
relation (see e.g. ICTBOSh . adopting the solar values and bolo- 
metric coiTections given in Table |8] For the primary a distance 
of 170.3 + 4.5 pc was found which corresponds well with the 
distance of 174.4 + 6.3 pc determined for the secondary compo- 
nent. For the system distance the average of both was adopted at 
172.0 + 3.8 pc. 



Notes. Bolometric corrections (BC) by [Flowed l ll996h were assumed 
together with T^q^ = 5780 K, BCq = -0.08 and M^^^ ^ = 4.74. 

Table 9. The effective temperature determined with five different 
methods. 



Method 


Primary 


Secondary 


Alonso calibration" 


6128d 


-60 


5205±86 


Holmberg calibration* 


6165d 


:68 


5248±109 


Casagrande calibration'' 


6244d 


:85 


5301±151 


Popper flux scale'' 


/ 




5098±58 


Spectroscopic temperature 


6100d 


:50 


5160±103 



References. 


Alonso etalj (Il996h 


ICasaerande et al 


(120 1 Oh '""iPoDoeild 19801) 



IHolmberg etall ( 120071) 



7. Stellar activity 

The light curves of EF Aqr show periodic variations that could 
be caused by star spots. The period of the variations differs, and 
the anomalies in the light curves disappear and reappear during 
the eight years over which the observations were made. Several 
subsets of the lightcurves were analysed separately, but it proved 
impossible to fit the variations. However, they are a clear indica- 
tor of the activity of one or both components. The Rossby num- 
bers, defined as the ratio of the rotation period over the convec- 
tive turnover time, for the primary and secondary components 
are approximately 0.71 and 0.08. The limit for higher magnetic 
ac tivity due to stronger dynamo action was determined at 0.65 
bv lHalll(ll994l) . This places the secondary well within the active 
region, and the primary at the border. The most compelling ev- 
idence for activity of both components is the Ca ii emission. In 
FigHthe Ca n - K (3933.68 A) and Ca ii - H (3968.49 A) emis- 
sion lines for both components are shown. In the spectrum of 
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Fig. 4. Ca ii-K (3933.68 A) and Ca n - H (3968.49 A) emission 
lines in the spectrum of both components. 

the primary, the emission lines appear in the hydrogen absorp- 
tion lines, while for the secondary, the Ca ii lines are visible as 
emission stronger than the continuum. These signs show that the 
primary is slightly active, while the secondary is very likely to 
exhibit strong chromospheric activity and star spots. 

Although the primary is less active than the secondary, it is 
prossible that starsports are the cause of the periodic variations in 
the radial velocity curve. Both components are suffuciently dif- 
ferent in temperature, mass and radius that the effect of starspots 
can differ strongly for both components. The primary being the 
larger and hotter component is likely to exhibit a stronger effect 
from starspots. 

8. Discussion 

In this section the absolute dimensions obtained for EF Aqr are 
compared with the properties of recent theoretical stellar evolu- 
tionary models. For an extended description of the models, the 
reader is refeiTed to CTB08. 

8.1. Comparison with solar-calibrated models 

Figures |5]-|7]illustrate the comparison with the Yonsei-Y ale (Y^) 
evolut ionary tracks and with the isochrones by Demara ue et al.l 
(I2OOI . The mixing length parameter in the convective envelopes 
was calibrated using the Sun, and was held fixed at l/Hp - 




6400 6200 6000 5800 5600 5400 5200 5000 



Teff [K] 

Fig. 5. EF Aqr compared with models for [Fe/H] = 0.00 (solid 
Une) and [Fe/H] = +0.10 (dotted Hne). Isochrones for 2.25, 3.5 
and 6.5 Gyr (left to right) are shown in black and tracks for the 
component masses are shown in red. 

1.7432. An enrichment law of Y = 0.23 + 2Z was a dopted, to- 
gether with the solar mixture of iGrevesse et alj (11996,) . leading 
to {X, YZ)q = (0.71564,0.26624,0.01812). Only models with- 
out a-element enrichment were considered ([a/Fe]=0.0). 

Fig. |5] shows that theoretical tracks for the observed masses 
and metalicity are hotter than the observations. Especially for 
the secondary, the temperature difference is large (~ 400 K). On 
the other hand, the theoretical temperature difference between 
the components of 800 K is only slightly smaller than the ob- 
served temperature difference of 965 + 125 K, which is partly 
covered by the uncertainty on the track positions caused by the 
small mass uncertainty. To incorporate the 0.10 dex uncertainty 
on the metalicity, theoretical tracks with [Fe/H] = -nO.lO were 
compared to the observed parameters, too. As was shown in Sec. 
|5] a higher metallicity gives rise to a higher spectroscopic tem- 
perature, leading to an increase of ~ 20 and ~ 40 K for the 
primary and secondary component. This higher metalicity track 
fits the observed temperature and radius of the primary within 
errors. For the secondary component, the theoretical parameters 
coiTespond better to the observed results, but the theoretical tem- 
perature is still 250 K hotter than observed, which is well beyond 
the uncertainty. 

The masses and radii calculated from the light and radial 
velocity curves are independent of the distance to the system 
(scale-independent). The mass-radius plane plotted in Fig. |6] 
shows that the model predicts an age difference of 4 Gyr be- 
tween the components, 2.25 Gyr for the primary and 6.5 Gyr for 
the secondary. Changing to isochrones with a higher metallicity 
makes no significant difference in this result. At the age of the 
primary, the predicted radius of the secondary would be 0.869 
Ro, which is 9 % smaller than observed, while the observations 
have a precision of The mass-luminosity relation, shown 
in Fig |7l predicts an age of around 2 Gyr for both components, 
with a metallicity of [Fe/H]=-i-0.10, or a younger age of around 
1 Gyr with a metallicity of [Fe/H] =0.0. 
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M/Mq 

Fig. 6. EF Aqr compared with models for [Fe/H] - 0.00 (solid 
line) and [Fe/H] = +0.10 (dotted line). Isochrones for 2.25, 3.5 
and 6.5 Gyr are shown. 



The same picture is seen i n the Victoria-Regina 
dVandenBerg et al.l l2006h and BaSTI dPietrinferni et all |2004 
models, which differ only slightly from the models with 
respect to the input physics and He enrichment laws. A short 
description of these models is given injCTBOS. 

In conclusion, the secondary component of EF Aqr is signif- 
icantly larger, and about 400 K cooler than predicted by mod- 
els, which adopt a mixing length parameter matching that of 
the Sun. The primary component agrees better with the mod- 
els. Both components, however, reveal a serious discrepancy be- 
tween models and the observations, which we attribute to prob- 
lems in the modelling, since the discrepancies are much larger 
than the observational uncertainty, which is seen particularly 
clearly in Figs|5]and|6] 

8.2. Comparison with mixing-iength tuned models 

Several authors have demonstrated that models that adopt a re- 
duced mixing length parameter are better at fitting active solar- 
type binary components. Both components of EF Aqr show clear 
signs of increased magnetic activity (in particular the secondary 
component), and have re latively high rotational velocities (> 10 
km s~ see e.g. ICBC09I) . making EFAqr a good candidate to 
compare to mixing-length tuned models. 

For this purpose, dedicated mod els for EE Aq r were com- 
puted with the GRANADA code by IClarell (120041) , which as- 
sumes an enrichment law of Y = 0.24 + 2.3Z t ogether with 
the solar mixture by Gre vesse & Sauval (Il998h . leading to 
(X,Y,Z)0=(O.7O4,O.279,O.O17). In this model, the mixing length 
parameter needed to reproduce the Sun is I /Hp - 1 .68. 

Several tracks for the precise observed masses were calcu- 
lated. Tests on the influence of the helium content were per- 
formed, and the helium content stated above gave the best re- 
sults. The five best tracks (three for the primary and two for 
the secondary) are listed with there resulting ages in Table [TOl 




M/Mq 

Fig. 7. EF Aqr compared with Y^ models for [Fe/H] = 0.00 (solid 
line) and [Fe/H] - -I-O.IO (dotted line). Isochrones for 2.25, 3.5 
and 6.5 Gyr are shown. 

Table 10. Ages of the primary and secondary components of 
EF Aqr determined from a comparison between the observed 
masses and radii with the Granada models for different mixing 
length parameters. 





Type 


l/H„ 


Age 


1 


P 


1.30 


1.50 + 0.18 


2 


P 


1.20 


0.90 + 0.20 


3 


P 


1.15 


0.70 + 0.20 


4 


S 


1.05 


1.55+0.60 


5 


S 


1.10 


2.25 ± 0.75 



The temperature-radius comparison between EF Aqr and the 
Granada models is shown in Fig. |8] The primary component is 
best fitted with the l/H,, -1.15 track at an age of 0.70 + 0.20 
Gyr, while the best fit for the secondary component is obtained 
with an l/Hj, - 1.05 track with resulting age of 1.55 + 0.60 
Gyr. When a slightly longer mixing length of l/Hj, - 1.30 is 
used, the age of the primary can be fitted at 1.50 + 0.18 Gyr. The 
correspondence with the temperature is slightly poorer, but still 
well within the uncertainty. Granada models with a solar mixing 
length parameter are unable to fit both components. 

The best-fitting tracks with identical derived ages for the two 
components are those with mixing length l/Hp - 1.30 for the 
primary, and I /Hp = 1.05 for the secondary. These tracks pre- 
dict a common age of 1 .5 + 0.6 Gyr for the system. This does 
imply that the mixing length parameter for both components can 
be tuned with reasonable precision, provided that accurate di- 
mensions and abundances are available. 

8.3. Comparison witli other binaries 

An overview of the properties of 1 1 solar-type eclipsing bina- 
ries of which the masses and rad ii are de termined with a pre- 
cision of at least 2% is given in ICBC09I . This list can be ex- 
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Fig. 8. EF Aqr compared with Granada models for the observed 
masses calculated with different mixing length parameters (see 
Table [TOl i. Primary component (filled cicles, upper left part of 
diagram): Case 1: full line, case 2: dashed line, case 3: dotted 
line. Secondary component (open circles, lower right part of di- 
agram): case 4: full line and case 5: dashed line. The circles on 
the evolutionary tracks mark ages of 0.5, 1, 1.5, 2, and 2.5 Gyr. 

tended with EW Orionis ( Clausen et all 1201 Oh and IM Virginis 
([Morales et al.l 12009.) . Clausen divided these binaries into two 
groups, the first consisting of active binaries with high rotational 
velocities, and the second group with longer periods and no signs 
of activity. (IM Vir belongs to the first group, while EW Ori be- 
longs to the second). The non-active binaries with a longer pe- 
riod can be fitted by solar-calibrated models, while the active 
binaries can be fitted with mixing-length-tuned models. EF Aqr 
is then a clear member of this second group, and supports the 
suggestion of Clausen et al. (2010) that the mixing length pa- 
rameter is governed by two main effects: For inactive main se- 
quence stars a slight decrease in mixing length with increasing 
temperature and mass, and secondly, for active stars a strong de- 
crease in mixing length compared with inactive stars of the same 
mass. 



9. Summary and conclusions 

From state-of-the-art observations and analyses, precise (0.6- 
1.0%) absolute dimensions have been obtained for the totally 
eclipsing GO system EF Aqr. A detailed spectroscopic analysis 
yielded a Fe abundance relative to the Sun of [Fe/H]=0.00+0. 10, 
with the abundances of Si, Ca, Sc, Ti, V, Cr, and Co and Ni being 
close to solar, too. 

The 0.956 M© secondary is found to be active with star spots 
and strong Ca ii H and K emission. The 1 .244 M© primary shows 
signs of activity as well, but at a lower level. The system orbit 
is circular, and the measured rotational velocities, 21+4 km s"' 
(primary) and 18+4 km s"' (secondary) excellently agreement 
with theoretically predicted rotational velocities for synchronous 
rotation. 



We showed that theoretical stellar models with solar- value 
mixing length parameters (Yonsei-Yale, Victoria-Regina and 
BaSTl) are unable to fit both components of EF Aqr at a com- 
mon age. At the derived age of the primary, the secondary was 
found to be 9% larger than predicted. Furthermore, the compo- 
nents were observed to be about 200 K (primary) and 400 K 
(secondary), cooler than predicted by the models. 

Models that adopt a significantly lower mixing length param- 
eter can remove these discrepancies. Granada models calculated 
for the observed metalicity with a mixing length of IjHp - 1.30 
(primary) and II Hp - 1.05 (secondary) reproduce the properties 
of both components of EF Aqr at a common age of 1.5+0.6 Gyr. 

EF Aqr falls in the group of solar-type systems with high 
rotational velocities that cannot be matc hed by solar-calibrated 
theoretical models as defined by'CBC09'. It supports the sugges- 
tion that magnetic activity is the likely cause for the radius and 
temperature discrepancies, which means in turn that active bina- 
ries have a significantly lower mixing length than their inactive 
counterparts of similar mass. 

This study is part of a larger project on solar-type eclipsing 
binaries: see CHOOl 
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